Subcritical dynamo bifurcation in the Taylor Green flow 
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We report direct numerical simulations of dynamo generation for flow generated using a Taylor- 
Green forcing. We find that the bifurcation is subcritical, and show its bifurcation diagram. We 
connect the associated hysteretic behavior with hydrodynamics changes induced by the action of 
the Lorentz force. We show the geometry of the dynamo magnetic field and discuss how the dynamo 
transition can be induced when an external field is applied to the flow. 
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Larmor [l[ is generally credited for suggesting that the 
magnetic field of the Sun (and, by extension, that of 
planets and other celestial bodies) could be due to dy- 
namo action - i.e. self-generation from the motions of an 
electrically conducting fluid. This principle has received 
much theoretical support 0] since then and has recently 
been validated by experimental observations 0, 0, H, S] . 
Dynamo action results from an instability: when the flow 
magnetic Reynolds number Rm exceeds a critical value 
R c Ml the null magnetic field state looses its stability to a 
non-zero magnetic field state. Because of the low value 
the magnetic Prandtl number of the considered fluids, 
this instability usually happens on a turbulent (noisy) 
basic state and the choice of an order parameter can be 
ambiguous Q- However, we can assume that the usual 
concepts of stability theory apply (cf. later) and study if 
the transition is supercritical or subcritical In most 
models and in all experiments, this bifurcation is super- 
critical: R C M is a unique number, albeit flow dependant. 



For instance R 



M 



14 and R C M 



18 for the constrained 



Karlsruhe and Riga experiments, while R C M ~ 32 for the 
fully turbulent VKS dynamo [5(. On the other hand, 
the dynamo bifurcation may also be subcritical partic- 
ularly because the action of a growing magnetic field is 
supposed to reduce hydrodynamic turbulence and main- 
tain dynamo action for lower Rm values. In fact, the 
transition can be globally subcritical if the basic state 
experiences instability with respect to finite amplitude 
perturbations A characteristic hysteretic behavior is 
then associated to the bifurcation, and the dynamo op- 
erates for a range of lower values R 9 M < Rm < R C M . 
Subcriticality has been discussed in MHD Alpha-Omega 
dynamical systems 1^, U | and also for numerical simula- 
tions of convective dynamos in spherical geometries [l2| 

In this Letter, we study the dynamo bifurcation using 
full MHD simulations, generated in a 3D-periodical do- 
main, by the Taylor Green forcing At low Reynolds 
numbers, this flow has several metastable hydrodynamics 
states [1J]. At higher Reynolds numbers, it has a well de- 



fined mean flow structure with superimposed intense tur- 
bulent fluctuations. Recent studies of the linear problem 
have shown that, while the dynamo thresholds may run- 



away in flows generated by random forcing [15j | , a dynamo 



is observed at all kinetic Reynolds numbers [l6|, [lj], EH 
in the Taylor-Green flow. We study the fully nonlinear 
regime and report here evidence of the subcriticality of 
the bifurcation. 

Using standard direct numerical simulation (DNS) pro- 
cedures wc integrate pscudospcctrally the MHD equa- 
tions in a 27r-periodic box: 



— + v • VB = B • Vv + 77V 2 B 

at 



-VP+jxB + fV 2 v + F, (1) 



(2) 



together with V • v = V • B = 0; a constant mass density 
p = 1 is assumed. Here, v stands for the velocity field, 
B the magnetic field (in units of Alfven velocity), j = 
( V x B)//io the current density, v the kinematic viscosity, 
r\ the magnetic diffusivity and V is the pressure. The 
forcing term F is given by the TG vortex 



' TG 



(M = 2/(t) 



sin(fco x) cos(fco y) cos(fco z) 
— cos(fco x) sin(fco y) cos(fco z) 




(3) 

implemented here at kg = 1. In the sequel, we consider 
two types of forcing: one in which f(t) is set to a constant 
- this is the constant force forcing (f(t) = 1.5) consid- 
ered in [ltij . In a second one 2/(t) is set by the condi- 
tion that the (1,1,1) Fourier components of the velocity 
remains equal to the Taylor-Green vortex - this is the 
constant velocity forcing considered in [13]. For the lin- 
ear instability problem, both forcing yield the same value 



of R C M [lfj, |17|. We now explore the non-linear regime, 



as well as the response to finite amplitude perturbations. 
Three control parameters drive the instability: the mag- 
netic and kinetic Reynolds numbers and the amplitude 



2 



of an external magnetic field Bo when applied. 
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(4) 



In the definition of the Reynolds numbers, the charac- 
teristic length scale is set to tt, the size of a TG cell 
when fco = 1. The characteristic speed v® ms is com- 
puted from hydrodynamic runs in which the Navier- 
Stokes equation is not coupled to the induction equa- 
tion, v® ma = {\j2Ev{t)) t - Here Ey is net kinetic energy 
Ey(t) and (•)( stands for averaging in time (1/T J T ■ dt). 
Likewise, in dynamo runs, the intensity of the magnetic 
field is estimated from the net magnetic energy EM(t), 
as b = (y/2E M (t)) t . 

Previous works [l?], EH have explored the response of 
TG flows to infinitesimal magnetic perturbations, as a 
function of the kinetic Reynolds number Ry. It was 
found that at any Ry, there exists a critical R C M above 
which perturbations grow exponentially This is illus- 
trated in Fig. [1] for a run at Ry = 563 and Rm = 281 
above the critical value R C M = 206. The initial magnetic 
field perturbation - with an energy level Em = 10 -17 - 
first grows exponentially. At time t ~ 300, the magnetic 
field has reached sufficient amplitude so that it can react 
back onto the velocity field, saturate the instability and 
reach a statistically stationary state, with approximate 
equipartition Em ~ Ey. Note that times are given here 
in units of equation (TT]), for which 1 is very close to one 
eddy turnover time of the flow (Tnl = T/ii, ms ~ 1.17). 
This transition from infinitesimal perturbations builds 
the (solid) red curve in Fig. 2. 

We have then quenched the system: at t = 1000, the 
magnetic diffusivity r\ is suddenly increased by a factor of 
4, lowering Rm below R M . After a short transient, both 
Ey and Em decrease and reach a second statistically sta- 
tionary state, with a non zero magnetic energy - a new 
dynamo state, for which equipartition is reached again 
(Fig. [T]). This behavior is an evidence for global sub- 
criticality The different levels of fluctuations in the 
two regimes suggest the possibility of different dynamo 
states, depending on the magnetic field or on history of 
the system. 

As subcritical bifurcations are also associated with hys- 
teresis cycles, we have repeated the quenching procedure 
starting from the same dynamo state A (obtained at 
t = 1000 at Ry = 563 in Fig. [I]) for increasing values of r], 
i.e. for decreasing Rm values. The (time-averaged) mag- 
netic and kinetic energy obtained after rearrangements 
are then recorded, and results summarized in Fig. [2] by 
the curve in the Bo = plane. Starting from point A, one 
can sustain the dynamo after quenching through points 
A2 to A9, until a value R 9 M substantially lower than R C M 
(at A9, Rm = 70 compared to Rm = 211 in A3). 

We have investigated further the system behaviour 
along the cycle by monitoring the spatial structure of 




FIG. 1: After a dynamo is self-generated from infinitesimal 
perturbations, the induction equation is quenched at t — 1000 
by a four-fold increase of the magnetic diffusivity. It corre- 
sponds to a sudden change from A to Ag - cf. Table I. 
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TABLE I: For each regime: root mean square a mplitude of 
the magnetic/ resp. velocity fields b = (y / 2Em(t)) , w rms = 
(y/2Ev(t)) , integral scale of the magnetic/resp. velocity fields 
L B = (T,E B (k,t))/k) , L v = (£Ev(k,t))/k) - E(k,t) is 
the uni- dimensional energy spectra. 



the magnetic and kinetic energies, so as to detect possi- 
ble changes in the flow structure. In a first regime, until 
point A7, the kinetic energy (and hence w rms , i.e. the 
turbulence intensity - see Table I) decreases and so does 
the magnetic energy - equipartition being essentially pre- 
served. Past A8, changes occur: Ey starts to increase 
abruptly, while Em continues to decrease, resulting in a 
decreasing ratio Em/Ev - see also Fig. 4. Other global 
quantities are also changing along this branch (see Table 
I). It corresponds to a modification in the spatial struc- 
ture of the magnetic energy. As can bee seen in Fig. [31 
the dynamo modes in A7 and ^48 arc different. At ^47, 
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the dynamo has a structure with magnetic energy 'tubes' 
in which the field line arc concentrated along diagonal di- 
rection (aligned with the energy structures). In A8, the 
dynamo has a magnetic energy with a wavy shape and 
the field line are no longer parallel to the energy struc- 
tures. In fact, the geometry of the A8 and A9 dynamo 
modes is reminiscent of the low kinematic mode of the 
TG dynamo [H. 



the kinematic low branch to the kinematic hig h branch 
as Ry increases (and turbulence develops) 17, 18]. The 
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FIG. 2: Bifurcation curves and hysteresis cycles when an ex- 
ternal magnetic field is applied (full diamond symbols) or 
without one (full circle symbols) . In this case, the subcritical 
quenched states (see text) form the red line. Jumps between 
the two branches link A to A' and C to C'. 





FIG. 3: Volume rendering (75% of max(6)) of the magnetic 
energy and magnetic field lines [l^J], for the normalized mag- 
netic field (B(pc,t)/B(t)) averaged in time during the run; 
(left) point A7 and (right) point A8. 

As turbulence influences the dynamo, we have re- 
peated the above sequence of quenching at varying ki- 
netic Reynolds numbers Rv- The result is shown in 
Fig. |H We first observe that the hysteretic behavior per- 
sists as Ry is lowered. In addition, the hysteresis cycle 
width, Rjy — R 9 M , decreases with Ry- It is interesting to 
compare their locations with respect to the dynamos win- 
dows evidenced in [l7|, EH for the Taylor-Green forcing. 
As shown in Fig.|4j R 9 M values are almost independent of 
Rv and lie close to the beginning of the first kinematic 
dynamo mode. Of course, the onset R C M switches from 



width of the dynamo cycle is thus linked to the evolution 
of the R c M (Rv) curve. 

The above results were obtained with a constant force 
scheme. We have also repeated the quenching procedure 
using the constant velocity forcing. As can be seen in 
Fig. 2] (black curve / diamonds symbols), the hysteretic 
behaviour remains, but the transition towards the non- 
dynamo state is more abrupt. Another difference con- 
cerns the response to quenching; with a constant velocity 
forcing we observed a lower magnetic saturation level b. 
Those differences could be explained by a change in hy- 
drodynamics properties such as the fluctuation level, at 
the same Reynolds number. In addition, when the veloc- 
ity is kept constant there may be less possibility for the 
Lorentz force to change the flow. 




FIG. 4: Hysteresis cycle for different Reynolds numbers and 
forcings - constant force (red, blue) and constant velocity 
(black). The thick solid line in the 6 = plane is the linear 
instability R C M vs Rv from dynamical runs; the kinematic 
dynamo windows [3], Rm £ [50,110] and Rai > 320, are 
delimited by the thick dotted lines. 

Finally, we have checked the influence of finite ampli- 
tude external perturbations on the hysteresis cycle by ap- 
plying an external magnetic field of amplitude Bq = 0.07 
in the vertical direction. The result at Rv = 563 is shown 
by the blue line in Fig. [2] When comparing to the Bq = 
case (red curve) , two effects are readily observed : (i) the 
hysteresis cycle is shortened and this is essentially due to 
a decrease in the onset R C M from infinitesimal perturba- 
tions; (ii) the amplitude of the magnetic energy in the dy- 
namo is decreased, as lower b values are obtained. These 
observations are indications that the external magnetic 
field has mediated a transition towards another equilib- 
rium state 1J] . The transition towards this second equi- 
librium state is quite robust: one can also obtain it by 
switching on the vertical magnetic field starting from a 
state with a well-developed dynamo (jump from A to A' 
in Fig. [5]). Conversely, starting from a dynamo state with 
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an applied magnetic field and switching it off, one returns 
to the zero-magnetic field hysteresis curve (jump from C 
to C' in Fig. 12). 

A less deterministic behaviour is observed when the 
system is operated in the vicinity of point D - shown 
along the blue curve in Fig. 2. At this point, the sys- 
tem is operated at a magnetic Reynolds number slightly 
smaller than the linear threshold (93.8 compared to 
about 100) and one observes that the the systems sponta- 
neously switches between dynamo and non-dynamo pe- 
riods, as shown in Fig [5] This is reminiscent of the "on- 
off" bifurcation scenario sometimes proposed for the dy- 
23| at high Ry- It has been observed in 
and experimental 12 51 1 versions of the Bullard 



namo 
models 124 

dynamo [2a] , and possibly in turbulent fluid dynamos Q . 
We note in Figj5]that the kinetic energy has stronger fluc- 
tuations during the dynamo periods. 




For the TG flow, we observe a decrease of the dynamo 
threshold by as much as 57 percent, with an external 
applied field of Bo = 0.07. We have also found that 
changes in the geometry of the dynamo states in the 
subcritical branch are consistent with the coexistence of 
several metastable hydrodynamics states (TJ]. Prelimi- 
nary observations in the VKS experiment also point to 
the existence of subcritical dynamos in the presence of 
global rotation j2?| , a feature also noted in some numer- 
ical models of the geodynamo [IH . 
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FIG. 5: Evolution on time of the kinetic (Ey) and magnetic 
energy (Eb) when the flow is operated in the immediate vicin- 
ity of point D - see Fig. 2. 

To summarize, we have evidenced in the TG flow sev- 
eral features characteristic of subcriticality of the dynamo 
instability. At variance with usual dynamical system, 
this behaviour is obtained in a fully turbulent system, 
where fluctuations are of the same order of magnitude as 
the mean flow. We may remark that in this case, the tra- 
ditional concept of amplitude equation may be ill-defined 
and one may have to generalize the notion of 'subcritical 
transition' for turbulent flows. Another feature is the sen- 
sitivity to perturbations of the order parameter through 
the application of an external magnetic field. The pertur- 
bation mainly acts through macroscopic changes in the 
system configuration (perturbation of the velocity field), 
allowing lower thresholds for dynamo instability. These 
findings open new perspective for experimental dynamos. 
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